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ANALYSIS OF BIDIRECTIONAL BALLOT SEQUENCES AND 
RANDOM WALKS ENDING IN THEIR MAXIMUM 

BENJAMIN HACKL, CLEMENS HEUBERGER, HELMUT PRODINGER, 

AND STEPHAN WAGNER 


Abstract. Consider non-negative lattice paths ending at their maximum height, 
which will be called admissible paths. We show that the probability for a lat¬ 
tice path to be admissible is related to the Chebyshev polynomials of the first or 
second kind, depending on whether the lattice path is defined with a reflective bar¬ 
rier or not. Parameters like the number of admissible paths with given length or 
the expected height are analyzed asymptotically. Additionally, we use a bijection 
between admissible random walks and special binary sequences to prove a recent 
conjecture by Zhao on ballot sequences. 


1. Introduction 

Lattice paths as well as their stochastic incarnation—random walks—are inter¬ 
esting and classical objects of study. Several authors have investigated a variety of 
parameters related to lattice paths. For example, Banderier and Flajolet gave an 
asymptotic analysis of the number of special lattice paths with fixed length in [2]. 
De Brujin, Knuth, and Rice |1] analyzed the expected height of certain lattice paths, 
and Fanny and Prodinger [T3] determined the asymptotic behavior of such paths 
with respect to several notions of height. 

The particular class of lattice paths we want to analyze in this paper is dehned as 
follows. 

Definition 1.1 (Admissible random walks and lattice paths). Let {Sk)o<k<n be a 
simple symmetric random walk on Uq or Z of length n starting at 0. That is, we 
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have P(S'o = 0) = 1 as well as 

^{Sk = j - 1 I Sk-i = j) = f‘{Sk = j + 1 I Sk -1 = j) = ^ for j > 1, 

F{Sk = 1 I Sk-i = 0) = 1, 

for random walks defined on Nq, and 

P(*S'fc = j - 1 I 5^-1 = j) = P(S'fc = j + 1 I Sk-i = j) = ^ for j e Z 

for random walks on Z. Then {Sk)o<k<n is said to be admissible of height h, if the 
random walk stays within the interval [0, h] and ends in h, i.e. Sk G [0, h] for all 
0 < fc < n and F{Sn = h) = 1. It is called admissible, if it is admissible of any height 
heN. 

The probability that a random walk of length n is admissible of height h is written 
as and for random walks on Nq and Z, respectively. Furthermore, the 
probabilities that a random walk is admissible at all are defined as pn '■= '^h>oP^'^ 
and Qn := respectively. 

Finally, an admissible lattice path is a sequence of integers realizing an admissible 
random walk. 

In a nutshell, this means that an admissible random walk is a non-negative random 
walk ending in its maximum. The definition is also visualized in Figure where all 
admissible lattice paths of length 5 are depicted. There are three admissible lattice 
paths of height 3, and one of height 1 and 5, respectively. Note that when considering 
random walks on Z, every lattice path has the same probability 2“”. Admissible 
random walks on Z are enumerated by sequence A167510 in m- 



Figure 1. Admissible lattice paths of length 5 

However, in the case of random walks on Nq, the probability depends on the num¬ 
ber of visits to 0: if there are v such visits (including the initial state), then the path 
occurs with probability 2“”+’'. Note that by “folding down” (i.e., reflecting about 
the a:-axis) some sections between consecutive visits to 0, or the section between the 
last visit and the end, 2'" lattice paths on Z can be formed, where the random walk 
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is never farther away from the start than at the end. We will call such lattice paths 
extremal lattice paths —and by construction, the number of extremal lattice paths of 
length n is given by To illustrate this idea of extremal lattice paths, all paths 

of this form of length 3 are given in Figure 



Figure 2. Extremal lattice paths of length 3 

One of our motivations for investigating admissible random walks originates from 
a conjecture in [16]. There, Zhao introduced the notion of a bidirectional ballot 
sequence: 

Definition 1.2 f [TUl Dehnition 3.1]). A 0-1 sequence is called a bidirectional ballot 
sequence if every prehx and suffix contains strictly more I’s than O’s. The number 
of bidirectional ballot sequences of length n is denoted by 3^. 

Bidirectional ballot sequences are strongly related to admissible random walks on 
Z. In fact, every bidirectional ballot sequence of length n-|- 2 bijectively corresponds 
to an admissible random walk of length n on Z: given an admissible random walk, 
every up-step corresponds to a 1, and down-steps correspond to 0. Adding a 1 both 
at the beginning and at the end of the constructed string gives a bidirectional ballot 
sequence of length n + 2. 

Therefore, bidirectional ballot walks may also be seen as lattice paths with unique 
minimum and maximum. 

While we restrict ourselves to simple lattice paths (i.e. the path has steps ±1), 
Bousquet-Melou and Ponty introduce a more general class of so-called culminating 
paths in |3|. Akin to bidirectional ballot walks, culminating paths are lattice paths 
with unique minimum and maximum—however, the lattice path has steps a and —b 
for hxed a, b > 0. They show that the behavior of these paths strongly depends on 
the drift a — b. In particular, for a = b = 1 (i.e. for bidirectional ballot walks) they 
determine the main term of the asymptotic expansion of Bn (cf. (3] Proposition 4.1]). 

In [16], Zhao also shows that Bn = Q(2^/n), states (without detailed proof) that 
Bn ~ 2”/(4n) and conjectures that 

Bn 1 1 

2”' 4n 6?7,2 \n^J 
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In this paper, we want to give a detailed analysis of the asymptotic behavior of 
admissible random walks. By exploiting the bijection between admissible random 
walks and bidirectional ballot sequences, we also prove a stronger version of Zhao’s 
conjecture. 

In order to do so, we use a connection b etwe en Chebyshev po lyno mials and the 
probabilities and (cf. Proposition 2.l| and Proposition 2.2, respectively), 
which we explore in detail in Section This allows us to determine explicit repre¬ 
sentations of the probabilities and qn-, which are given in Theorem [Tj The analysis 
of the asymptotic behavior of admissible random walks of given length shall focus in 
particular on the height of these random walks. In this context, we dehne random 
variables Hn and Hn by 


nHn = h) : = 


(h) (h) 

15 ^. P(A = h) := 

Pn qn 


These random variables model the height of admissible random walks on No and Z, 
respectively. Besides an asymptotic expansion for and qn, we are also interested 
in the behavior of the expected height and its variance. The asymptotic analysis of 
these expressions, which is based on an approach featuring the Mellin transform, is 
carried out in Section and Section and the results are given in Theorem and 
Theorem]^ respectively. Finally, Zhao’s conjecture is proved in Corollary |4.2[ 


2. Chebyshev Polynomials and Random Walks 

We denote the Chebyshev polynomials of the hrst and second kind by and Uh, 
respectively, i.e., 

Th+i{x) = 2xTh{x) - Th-i{x) for h > 1, Tq{x) = 1, Ti{x) = x, 

Uh+i{x) = 2xUh{x)-Uh-i{x) for h > 1, Uo{x) = l, Ui{x)=2x. 

In the following propositions, we show that these polynomials occur when analyz¬ 
ing admissible random walks. As usual, the notation [z^]f{z) denotes the coefficient 
of z^ in the series expansion of f{z). 

Proposition 2.1. The probability that a simple symmetric random walk {Sk)o<k<n 
of length n on Ij is admissible of height h is 

(1) = P(0 <So,Si,...,Sn<h and Sn = h) = 2[z^+^]—^j— 

for h > 0 and n > 0. 






BALLOT SEQUENCES AND RANDOM WALKS 


5 


Proof. We consider the (h + 1) x {h + 1) transfer matrix 


/O 

1 

2 


M, = 


1 

2 

0 


0 

1 

2 


0^0 


\0 0 


0\ 

0 

0 


1 

2 

0 / 


which has the following simple yet useful property: if wl^l is the probability that 
0 < Sq, Si,..., Sn < h and Sn = k, then the following recursion for the vectors 
= {Wnl ' ' ' ’ hoHs: 


■ Mh = w. 


(h) 
n+l 1 


SO ■ MJ^. The initial vector is = Cq = (1,0,... ,0). Since we also 

want that Sn = h, we multiply by the vector e/i = (0,..., 0,1)"'^ at the end to extract 
only the last entry This yields the generating function 


P) _ ^ _ 


nW’yU — 
rC 

Cramer’s rule yields 


= ^eoMfchZ^ = eo(/ 


zMh) 




hn — h 


z'^2 


n>0 


det(J — zMh) 


The determinant of / — zMh can be computed recursively in h by means of row 
expansion, see (for instance) [U p.97]: 

det(/ - zMh+ 2 ) = det(/ - zMh+i) - ^ det(/ - zMh). 

Comparing this with the recursion for the Chebyshev polynomials and checking the 
initial values, we hnd that - —Therefore, we obtain 




{h)^n _ 


n>0 


zUh+iil/z) 

from which ([^ follows by extracting the coefficient of 


□ 


An analogous statement holds for admissible random walks on Nq with the sole 
difference that in this case, the Chebyshev polynomials of the hrst kind occur. 
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Proposition 2.2. The probability that a random walk {Sk)Q<k<n of length n on No 
is admissible of height h is given by 

(2) = P(0 <So,S^,...,Sr,<h and Sn = h) = ^ 


Th+iO-/z) 


for h>D and n > 1. 


Proof. For random walks with a reflective barrier at 0, the (h + 1) x (h + 1) transfer 
matrix has the form 


M, = 


/o 

1 

0 

. 0\ 

1 

2 

0 

1 

2 

. 0 

0 

1 

2 

0 

0 




1 




, ■ 2 

\0 

0 

0 

... 1 0 ; 


By the same approach involving Cramer’s rule as in the proof of Proposition 2.1 
And the generating function 

^h2i-h 


we 


J2Pn^^" = ^0{I-zMk)-W = 


n>0 


det(/ — zMh) 


where we have the recursion 


det(/ - zMh+ 2 ) = det(/ - zMh+i) - — det(/ - zMh) 
for the determinant of / — zMh- Finally, (2) follows from -— 

which can be proved again by verifying that the same recursion holds for the Cheby- 
shev-T polynomials and that the initial values agree. □ 


Remark. The coefficients of 


Th{ijz) also been studied in [9]. There, the case 

of fixed h is investigated, whereas we mostly focus on the asymptotic behavior of 
n ^ 00. 


Using the results from Proposition |2.1 and Proposition 2.2, we may give explicit 
representations of the probabilities and by investigating the Chebyshev poly¬ 
nomials thoroughly. 

Remark (Iverson’s notation). We use the Iversonian notation 

1 if expr is true. 


|expr] = 

popularized in P Chapter 2]. 


0 otherwise. 
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In the following theorem and throughout the rest of the paper, m will denote a 
half-integer, i.e., m G = {|, 1, |, 2,..While this convention may seem unusual, 
it simplihes many of our formulas and is therefore convenient for calculations. 


Theorem 1. With Th,k '■= (h + l){2k + l)/2 and Vh,k '■= {h + 2){2k + l)/2, we have 


( 3 ) 

{h) 

P2rn-1 = 

( 4 ) 

Qp) - 

H2m-2 — 


E(-i) 


fc >0 


kWk 
m \m 


4 ‘^vl k-m 


4™ ( 2 m — l)m \m — Vh,k 


2m 

Th,k 

2m 


|h -I- 1 = 2m mod 2 ], 
|h = 2m mod 2 ] 


for h > 0 and half-integers m G with m > 1. 

Proof. We begin with the analysis of The probabilities are related to the 

Chebyshev-T polynomials by Proposition |2.2[ It is a well-known fact (cf. [I2l 22:3:3]) 
that these polynomials have the explicit representation 


Th{x) = 

which immediately yields 

( 5 ) 


{x — \/x^ — 1 4 - (x - 1 - — 1 


Thil/z) 


{1 - + {1 + 


=: z’^Y{z‘^). 


By applying Cauchy’s integral formula, we obtain the coefficients of the factor Y[t) 
encountered in (|^. We choose a sufficiently small circle around 0 as the integration 
contour 7 . Thus, we get 

[r]F(t) = ^ 




1 

27ri 


(1-Vi-0^ + (i + \/i-0^ 


dt. 


We want to simplify the expression this integral. This can be achieved by 

the substitution t = -( 77 ^, which gives us dt = {1 — u) ■ 7177 ^ du and a/I — t = 

Also, the new integration contour is 7 , which is still a contour that winds around 
the origin once. Then, again by Cauchy’s integral formula, we obtain 

+ 1 


inn<) = ^ui-«) 


= |'U "](1 - u) 


22n+h-l('t ^n+l 

(1 + 


du 


22 n+h-l^l J^uhy 
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Expanding the factor — into a series with the help of the geometric series 

and the binomial theorem yields 


(1 + n) 


2n-\-h—l 


2n-\-h—l 


1 + m'’ 


and therefore 


^(- 1 )= ^(- 1 )^ 

k>0 k>0 j=0 


2n + h — 1 

j 




m 


(l + n) 


2n-\-h—l 


l + u^ 




fc >0 


2?7, + h — 1 
i-hk 


This allows ns to expand the expression enconntered before, that is 


[i”]r(t) 


[“” 1(1 - 
1 

22n+/i—1 


+ 

22n+h-l(^l 



/2n + h — 1\ 
\n — hk — 1J 


Using the binomial identity 

'N-1 
a 


(N -1\ _ N -2a /iV\ 
Va- ly N Vay’ 


the expression above can be simplihed so that, together with ([^, we hnd 
1 


n(i/z] 


= 2 


f ^ \ 2 n+h 2hk -\- h f 2n + h 

2^ \ 9 / 2-^ 2n + h \n — hk 


n>0 


k>0 


By plngging this into ([^, we obtain 

p('>> = 2 [.-“+‘] ^ 


Th+ii^/z) 

z^2i+h+i^^ ^^^2{h + l)k + h + l f 2i + h + l 


=4[--""iE(D E(-i) 


2 / 1-1 


l>0 

„n—/i-j ^ ^ 

.2 


k>0 


2i + h + l \i-{h + l)k 




;-2(h + l)fc + h + l / 2i + h + l 


e>o 


k>0 


2i + h + l \i-{h + l)k 


Combinatorially, it is clear that = 0 for n and h of different parity, as only 
heights of the same parity as the length can be reached by a random walk starting 
at the origin. This can also be observed in the representation above. Assnming 
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n = h mod 2, we can write n — h = 2i or equivalently = i. This gives us 


^ 1 + 1 )^ + ^ + 1 
^n-1 2-J^ ^ 


k>0 


n + 1 




fc >0 


n + 1 -{h + l)k^ 

{h + l)(2/c + 1) / n + 1 

^Ti - \{h + i){2k + i))' 


Substituting n = 2m — 1 with a half-integer m G such that h + 1 = 2m mod 2, 
and recalling that Th^k = {h + l){2k + l)/2, the representation in ([^ is proved. 

For the second part, we consider the explicit representation 


Uh{x) = 


{x + — 1 — (x — \/x2 — 1 


2 ^ 


x^ 


of the Chebyshev-f/ polynomials, which is equivalent to 


= z 


2 ^/^ 


Uh{)-/z) (1 + — (1 — \/l — )^+i 

Formula Q is now obtained in the same way as ([^. 


□ 


With explicit formulae for the probabilities and q^\ we can start to work 
towards the analysis of the asymptotic behavior of admissible random walks. 


3. Admissible Random Walks on Nq 

In this section, we begin to develop the tools required for a precise analysis of the 
asymptotic behavior of admissible random walks on Nq- 

Recalling the result of Theorem we find that in the half-integer representation 
P 2 m-i) fhe shifted central binomial coefficient appears. Hence, for the pur¬ 
pose of obtaining an expansion for P 2 m-i = analyzing the asymptotics 

of binomial coefficients in the central region is necessary. In the following, we will 
work a lot with asymptotic expansions. The notation 

OO 

e=-L 

(as n —)■ cx)) is understood to mean 

R-l 

f{n) = ^ aen~^ + 0{n~^) 

£=-L 
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for all integers R > —L, even if the series does not converge. Likewise, an asymptotic 
expansion in two variables given by 


oo J(£) j 

f{a,n) ~ 

l=-Lj=0 

is to be nnderstood as 

R-i JW j 

/(a,n) = + 

e=-Lj=o ^ 

for all R > —L. 


Lemma 3.1. For n G and |a| < such that n — a G N, we have 


2n 

n — a 

2j 


rni 


a 


exp (- - S{a, n) 


n 


with S{a,n) ;= and 


(6) c,j = [a^^n ']( E 


dr 


r>0 


r>0 


i-iydr 

{n + ay 


(-lyd, 


Y — 

(r) — 


X 




r>0 


-l/2\ a 


r In 


.,2r 


T 2 r 


r>0 

1 


{n — ay 


E l lx ^—\ 

r^ n^^ ^ 


r>0 


a 


2i 


L t>0 


{t + 2){2t + 3) n‘^\ 


where the coefficients dr come from the higher-order Stirling approximation of the 
factorial, cf. Q. Additionally, the estimate 

S{a,n) = l + o(^^^ 

\ n 


holds for |a| < and we know that cqo = 1 as well as cij = 0 if j > 
If \a\ > nf!'^, the term 

( V 4 " = 0(exp(-n^/3)) 

\n — a J 

decays faster than any power of n. 


Proof. We begin by recalling the higher-order Stirling approximation (cf. [71 p. 760]) 


( 7 ) 
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An explicit representation of the coefficients dj can be found in [10]. From the 
logarithmic representation of the factorial (see [T] p. 766]), the expansion 


( 8 ) 


n\ A/2vm Vn 


E 

j>0 




for the reciprocal factorial follows. 

Let us assume |q!| < Then, by applying ([^ and ([^ to the shifted central 

binomial coefficient, we obtain 

2n \ {2n)\ 


n — aj {n — a)\{n + a)\ 
1 


2n 


n? J (n + a)"'+"(n — aY 


(2n) 


X 




5:(lu)(5: 


r>0 


r>0 


i-iydr 

(n + aY 


^ {n — a) 


r>0 


The factor (l — can be expanded as a binomial series, resulting in 


1 -^ 


q,2x -1/2 






r>0 


-l/2\ a 


r n 




,2r 


The remaining factor is handled by means of the identity n” = exp(nlogn), which 
leads to 

(2n)^" 

-—;———-r— = exp(2nlog(2n) - {n + a) log(n + a) - {n - a) log(n - a)) 

(n+ (n — a)" “ 

= exp(2n log 2 + 2n logn — (n + Q!)(logn + log(l + a/n)) 

— {n — a)(logn + log(l — a/n))) 

= d"' exp(a log(l — a/n) — a log(l + a/n) 

— nlog(l — a/n) — nlog(l + a/n)). 

By expanding the logarithm into a power series, we can simplify this expression to 


(2n) 


2n 


(n + Q;)”+"(n — aY 


= 4”exp 2 


-E 


1 a 


2i+2 




1 a 


21+2 


2f + 1 ^ 2f + 2 71^*+^ 

t>o t>o 

/ a^\ / a^ 1 a^* 

= 4 .xp(--)exp( 


a 
n 

0? 

= 4^ exp (- 

n 


1 r 


r>0 


a 


2t 


t>0 


(t + 2)(2t + 3) n^Y 





















12 


B. HACKL, C. HEUBERGER, H. PRODINGER, AND S. WAGNER 


We also use 

1 11 1 / Cl \ ^ 

n±a 77,1 it - nV nJ 

n r>0 

By the symmetry of the binomial coefficient, the resulting asymptotic expansion has 
to be symmetric in a. Assembling all these expansions yields the asymptotic formula 



where S{a,n) is dehned as in the statement of the lemma. 

Note that do = 1) cind thus the hrst summand of the series in ([^ is 1—which gives 
Coo = 1- Summands where the exponent of a exceeds the exponent of 1/n only occur 
in the last series, with the maximal difference being induced by Thus, if 

j > 1^, we have C£j = 0. Together with |ci| < this implies the estimate for 
S{a, n). 

For |ci| > we can use the monotonicity of the binomial coefficient to obtain 

/ 2n \ / 2n \ 

\ji — a) ~ \?7, — / ’ 

for which the exponential factor ensures fast decay, 

exp (-h^)=0 (exp (-„'/»)), 

and as everything else is of polynomial growth, the statement of the lemma follows. 

□ 


Now that we have an asymptotic expansion for the shifted central binomial coef- 
hcient, let us look at our explicit formula in ([^ again: we have 


{h) 

P2rn-1 = 


E(-i) 


k 


A:>0 


2m 


m \m- Th,k 


■ |h + 1 = 2m mod 2], 


where Th^k = {h + l){2k + l)/2. Therefore, the total probability for a random walk 
of length 2m — 1 on Nq to be admissible is given by 


P2m-1 — 


E 

h>0 


(h) 

P2rn-1 = 


h,k>0 

h+l=2m mod 2 


2m 

m \m- Th,k 


The terms where can be neglected in view of the last statement i 


m 


Lemma 3.1, as their total contribution decays faster than any power of m: note that 
there are only 0{m?) such terms (trivially, h, k < m), each of which contributes 
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0(mexp(—to the sum. For all other values of h and k, we can replace the 
binomial coefficient by its asymptotic expansion. This gives us, for any L > 0, 


P2m—1 





h-\-l=2m mod 2 


+ 0 



E 

h-\-l=2m mod 2 


2J(L)+1 
‘ h,k 
^L+l 




L -1 


e=o j>o 


'^h,k 


where J{L) < |L since Qj = 0 for j > Since the sum clearly contains 
terms, the error is at most The expo¬ 

nent can be made arbitrarily small by choosing L accordingly. Finally, if we extend 
the sum to the full range (all integers h, k > 0 such that h -|- 1 = 2m mod 2) again, 
we only get another error term of order 0(exp(—which can be neglected. In 
summary, we have 


( 9 ) 


P2m—1 


E 


Jmii 

^ h,k>0 

h+l=2m mod 2 


2 2 j 

T^k'^h,k ( '^h,k\ '^h,k 

If—exp- -) 2^ 

m \ m / 

Ci>o 


This sum can be analyzed with the help of the Mellin transform and the converse 
mapping theorem (cf. [B]). In order to follow this approach, we will investigate those 
terms in ([^ whose growth is not obvious more precisely. That is, we will focus on 
the contribution of terms of the form 


E (-1)' 


-2j+l . 
‘ h,k 


h^k>0 

h-\-l=2m mod 2 


T-2 

T^h,k 


m 


We are also interested in the expected height and the corresponding variance and 
higher moments of admissible random walks. Asymptotic expansions for these can be 

(h) 

obtained by analyzing moments of the random variable Hn with P(i7n = h) := 
as stated in the introduction. For the sake of convenience, let us consider the r-th 
shifted moment E(if 2 m-i + 1)^- We know 

+ 1)’' = + l)’'P(i/2m-l = h) = _ 


The asymptotic behavior of the denominator is related to the behavior of the sum 
from above—and fortunately, the behavior of the numerator is related to the behavior 
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of the very similar sum 


E (-lt^^‘(A+irexp 

h,k>0 

h+l=2m mod 2 



The following lemma analyzes sums of this structure asymptotically. 


Lemma 3.2. Let j, r G No- Then we have 

(10) E ‘(A+f “p (- 

h^k>0 

/i+l=2m mod 2 

= + 1 + ^ 0{m-^) 

for any fixed K > 0, where (3{-) denotes the Dirichlet beta function. 


Remark. The Dirichlet beta function is also often called Catalan beta function, and 
it is dehned by 

^ f _1 \k 

'"'■SifeV 

It can be expressed in terms of the Hurwitz zeta function as (3{s) =4 ^(C(s, 1/4) — 
C(s,3/4)). Amongst many other interesting properties, it satishes the zeta-like func¬ 
tional equation (cf. |T2], 3:5:2]) 

ffil — s) = (7r/2)“^ sin(7rs/2)r(s)/9(s), 

which also implies that /d(s) has zeros at all negative odd integers. 


Proof of Lemma \3.S\ If we substitute m = x the left-hand side of (10) becomes 

/(!):= E + 

h^k>0 

h-\-l=2m mod 2 


This is a typical example of a harmonic sum, cf. [6l §3], and the Mellin transform 
can be applied to obtain its asymptotic behaviour. First of all, it is well-known that 
the Mellin transform of a harmonic sum of the form f{x) = Y2k>i ^kdib^x) can be 
factored as J2k>i^kb'^^g*{s) [6l Lemma 2], provided that the half-plane of absolute 
convergence of the Dirichlet series A(s) = non-empty intersection 

with the fundamental strip of the Mellin transform g* of the base function g. In this 
particular case, the Dirichlet series is 


A(s):= E + 

h,k>0 

h-\-l=2m mod 2 
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and the base function is g{x) = exp(—with Mellin transform g*{s) = |r(|) and 
fundamental strip ( 0 , oo). 

Now we simplify the Dirichlet series. For s G C with 3fJ(s) > 2j + 2 + r, the sum 
A(s) = 

/i,fc>0 

/i+l=2m mod 2 

converges absolutely because it is dominated by the zeta function. In view of the 
definition of the (3 function, this simplihes to 

A(s) = - (2j + l))K 2 m(s - (2j + 1 + r)), 


where K 2 m(s) depends on the parity of 2m. We hnd 


^2m(s) — 'y ^ (h + 1) 

h>0 

h+l=2m mod 2 


2 ®C('S) for m eN, 

(1 — 2 “^)C(s) for m ^N. 


Thus, the Mellin transform of / is 

r(s) = A(ii)j-(s) = ir(02*-P^+‘y(s - (2j + l))^ 2 „(s - (2j + 1 + r)). 

By the converse mapping theorem (see Theorem 4]), the asymptotic growth of 
f{x) for a: —)■ 0 can be found by considering the analytic continuation of f*{s) further 
to the left of the complex plane and investigating its poles. The theorem may be 
applied because A(s) has polynomial growth and r(s/2) decays exponentially along 
vertical lines of the complex plane. 

We find that f*{s) has a simple pole at s = 2j + 2 + r, which comes from the zeta 
function in the definition of K 2 m- There are no other poles: f3 is an entire function, 
and the poles of T cancel against the zeros of (3 (at all odd negative integers, see the 
earlier remark). 

The asymptotic contribution from the pole of f* is 
Res(r, s = 2j + 2 + r) • + 1 + 02"++ l)0-(2j'+2+") 

= 2 "-W(j + 1 + 0/3(r + l)m^+i+"/2, 


m 


which does not depend on the parity of 2 m, as the respective residue of A 2 m is 
either case. Finally, the O-term in (10) comes from the fact that f* may be continued 
analytically arbitrarily far to the left in the complex plane without encountering any 
additional poles. □ 


Remark. In Lemma 3.2, particular values of the Dirichlet beta function are required. 
To compute the asymptotic expansions for the first moments, we need /3(1) = vr/d. 
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/3(2) = G ~ 0.91597, as well as (3{3) = 7r^/32, where G is the Catalan constant. 
These values are taken from [121 Table 3.7.1]. 


At this point, all that remains to obtain asymptotic expansions is to multiply the 
contributions resulting from Lemma 3^ with the correct coefficients and contribu¬ 
tions from (|^. 

Theorem 2 (Asymptotic analysis of admissible random walks on No). The proba¬ 
bility that a random walk on Nq is admissible can be expressed asymptotically as 

TT bV^TT i2(\/27i ib(i\/2Tr i»ypyi^v:^vr 1 

2n 


( 11 ) pn = — 



^ 127^^ 1571\/^ 1896913^^ 




'n 


11 


24^^ 960\/^ 16128v/n^ 184320vV 

where \/'n'l2 ^ 1.25331. The expected height of admissible random walks is given by 


2n 

(12) E77„ = 2GW — 


5v/2G 13iy2G 1129^2^ 

1 + -^ ^ 


60 ™ 720V 


+ o( 


12096-0™^ 

where 2G-\/2/7r 1.46167, and the variance of Hn can be expressed as 


(13) YHn = 


TT^ - 32G2 TT^ - 40G2 
-n -\- 


471 


67r 


TT^ - 12G^ llvr^ - 265G2 ^ / 1 

+ —7WXX—1— + 0 — 


1807m 


18907™^ 




where (tt^ — 32G^)/(47r) 
given by 


(14) 


eh: 


0.33092. Generally, the r-th moment is asymptotically 
2 »’/ 2-|-2 


TT 


T(^ + l)/3(r + l) 02 _ 


Moreover, if rj = h/^/n satisfies 3/0ogn < rj < yJ\ogn/2 and h = n mod 2, we 
have the local limit theorem 


Ah) 

(15) f>{Hr. = h) = ^ 

P'^ 


(16) 


ph ^ 20(7]) 8r] 


n 


Tix/n 


5 ^(-l)'=( 2 /c + l)exp 


k>0 


{2k + l)W 


2^ 


rj" 


^(-l)^(2fc + l)exp 


k>0 


7l\2k + lf 
8rf 


Remark. The fact that the two series in (15) and (16) that represent the density 0(7^) 
are equal is a simple consequence of the Poisson sum formula. We also note that 
the asymptotic behavior of the moments of Hn readily implies that the normalized 
random variable Hn/^/n converges weakly to the distribution whose density is given 
by 0 ( 77 ) (see |3 Theorem C.2]). The local limit theorem (15) is somewhat stronger. 

Proof. With ([^ and the result of Lemma 3.2, obtaining an asymptotic expansion of 
P 2 m-i is only a question of developing the shifted central binomial coefficient and 
multiplying with the correct growth contributions from (10). By doing so (with the 
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help of SageMath [15]: a worksheet containing these computations as well as some 
numerical comparisons can be found at http://arxiv.org/src/1503.08790/anc/ 
random-walk_NN. ipynb), an asymptotic expansion in the half-integer m is obtained. 
Substituting m = (n -|- l)/2 then gives 0- 


The results in (12) and (13) are obtained by considering 


E{Hn + !)'■ = 


Efe>o(^ + 

Pn 


(h) 


making use of ([^ and Lemma [3.2| again. Note that we have EH^ = IE(i7„ -|- 1) — 1, 
as well as = E[Hn + 1)^ — [E{Hn -l- 1)]^. For higher moments, we only give the 
principal term of the asymptotics, which corresponds to the coefficient cqo in (|^, but 
in principle it would be possible to calculate further terms as well. 

It remains to prove (15). To this end, we revisit the explicit expression (recall that 
we set n = 2m — 1) 


(h) 

P2rn-1 


4 

4™ 




'^h,k 

m 



First of all, we can eliminate all k with since their total contribution is at 

most 0(mexp(—as before. For all other values of /c, we replace the binomial 


coefficient according to Lemma 3T by 

4 m 


m 


exp 


2m 

m - Th,k 


a/FtF 


exp 


!M)(i + o 

m 


1 + 'Th,k 


m 


Note here that 

and likewise 

It follows that 
Th,k 


'Th,k = 


{h + l){2k + l) h 


— — (2A;-|-l)|l-|-0( — 


h 


‘h,k 

m 


h‘^{2k + lf 
2n 




h,k 

m 


h{2k + l) 


n 


exp 


l + 0(w + - 
n 


h‘^{2k + lf 
2n 


i + o(i + lF±l 

ri n 


We are assuming that r^^k < m^^^ = ((n-|-l)/2)^/^, which implies /h). 

In view of our assumptions on h, this means that the error term is 0{n~^/^^J\ogn). 
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Thus we have 


Ah) _ (h) _ 

Pn ~ P2m-1 ~ 


A^/2h 




X 


k>0 

T'h,fe<((n+l)/2)2/3 


(-l)'=(2A: + l)exp (- 


h\2k + iy 

2n 


1 + 0 


A/logu 

uV 6 


0(nexp(-(n/2)^/3)). 


Adding all terms Th^k > = {{n + l)/2)^0 back only results in a negligible 

contribution that decays faster than any power of n again, but we need to be careful 
with the 0-term inside the sum, as we have to bound the accumulated error by the 
sum of the absolute values. We have 

h^{2k + lf' 


Y^(2k + 1)exp 


fc >0 


2n 


= 0{n/h^ 


which can be seen e.g. by approximating the sum by an integral (or by means of the 
Mellin transform again), so 




'7171'^ 


-l)\2k + l)exp 


k>0 


AV2ri 


Tin 


5^(-l)'^(2fc + l)exp(- 


k>0 


Since Pn = \/^(l + 0{n ^)), this yields 
JA) Q„ 

” 5^(-lf(2i+l)exp 


Pn 

Pr 


Ti\/n 


k>0 


h‘^{2k + lf 
2n 

h'^{2k + lf 
2n 


p^{2k + 1 )" 


+ o( 


+ 0 


hn^/^ 

logn 


n 


n 


7/6 


+ o( 


logn 


n 


2/3 


20 (h) 


+ 0 


n 


logn\ 


For h ^ 1) the sum is bounded below by a constant multiple of exp(—r/^/2) (as 
can be seen by bounding the sum of all terms with > 1), which in turn is at 
least exp(—(logn)/8) = by our assumptions on p. Thus the hrst term in¬ 

deed dominates the error term in this case. If h < 1, we use the alternative repre¬ 
sentation (16), which shows that 0(h) is bounded below by a constant multiple of 
h“^exp(—7r^/(8h)^). This in turn is at least (l/9)n“’^ /^^logn by the assumptions 
on p, and since 7r^/72 < 1/6, we can draw the same conclusion. □ 
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Remark. As stated in the introduction, the number gives the number of extremal 
lattice paths on Z—and thus, with the asymptotic expansion of Pn, we also have an 
asymptotic expansion for the number of extremal lattice paths on Z of given length. 

This concludes our analysis of admissible random walks on Nq. In the next section, 
we investigate admissible random walks on Z. 


4. Ballot Sequences and Admissible Random Walks on Z 

In principle, the approach we follow for the analysis of the asymptotic behavior of 
admissible random walks on Z is the same as in the previous section. However, due 
to the different structure of Q, some steps will need to be adapted. 

With the notation of Lemma 13.11 we are able to express q 2 m -2 for a half-integer 


m G with m > 1 as 


(17) 


q2m-2 






m 


rrmr 2m — 1 ' m 

h,k>0 

h=2m mod 2 


exp ( — 


V 


h^k 

m 




V 


2j 

h,k 


e,j>o 




In analogy to our investigation of admissible random walks on No, we also want to 
determine the expected height and variance of admissible random walks. These are 
related to the random variable which we dehned by 

_ (h) 

P(R„ = h) = ^. 

Qn 

To make things easier, we will investigate moments of the form E,{Hn + 2)^. They 
can be computed by 


E{Hn + 2Y = + 2YF{Hn = h) = 


h>0 


Xfe>o(^ + 

Qn 


(h) 

n 


Therefore, we are interested in the asymptotic contribution of 


E 


m 


h,k>0 

h=2m mod 2 


m 




which is discussed in the following lemma. 

Lemma 4.1. Let K > 0 be fixed. Then we have the asymptotic expansion 


(18) 


E 


- m 


h,k>0 

h=2m mod 2 


m 


exp ( — 


V 


h,k 

m 


'miT 


+ 0{m 
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For j G N we have 


(19) 


E 


27)2 


m 


h^k>0 

h=2m mod 2 


m 


^h,k exp ( - 


V 


h,k 

m 


/l0£772 1 / i\ 1 \ 

= + 27 + log2+ -ijjyi + -j + — + [m ■ (21og2 - 2)j 

X + 0{m~^) 

where iIj^s) is the digamma function. Finally, for j G No? t G N we find 


( 20 ) 


E 


27)2 

‘^^h.k 


m 


l'»i(fc + 2)'exp( - 


h,k>0 

h=2m mod 2 


m 


V 


h^k 

m 


= ir(j + 

r + 1' 


+ 2^3 + “ l)C(r’ + ^0(m ^), 

where K 2 m(s) = 2“*C('S) for m G N and ti 2 m{s) = (1 — 2“^)('(s) — 1 otherwise. 
Proof. Let j, r G Nq. We want to analyze the sum 


E 


27)2 


m 


vl’,(h + 2Ye^p(- 


h=2m mod 2 


m 


V 


h,k 

m 


asymptotically, where m is a half-integer in with m > 1. 

In analogy to the proof of Lemma 3.2, we substitute x~‘^ = m, so that the sum 
becomes 


fix) ■= Y1 + 2 )" 


h,k>0 

h=2m mod 2 


2^' E dir{h + 2reM-vUA- E I’h('" + 2)''exp(-t-7x2 


h,k>0 

h=2m mod 2 


h,k>0 

h=2m mod 2 


=: 2x^fi{x) - f2ix). 


Both /i and /2 are harmonic sums, and we determine their Mellin transforms as 
we did earlier in the proof of Lemma [3^ By elementary properties of the Mellin 
transform, we know that f*{s) = 2/^(s + 2) — (s). Let Ai and A 2 be the Dirichlet 
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series associated with the harmonic sums /i(x) and / 2 (t), respectively. We hnd 

Ai(s) = ^ + 2y = 2"-(2t+2) ^ (h + + l)2i+2-- 

h,k>0 h,k>0 

h=2m mod 2 h=2m mod 2 

= (2"-(2i+2) - l)C(s - (2j + 2)) + 

/i>0 

h=2m mod 2 


We investigate the sum over h separately, and obtain 

^2m(s) := (h + 2) * = 


/i>0 

h=2m mod 2 


2-‘C(s) 

(1-2-KM-i 


for m G N, 
for m yN. 


Therefore, we hnd the Mellin transform of the hrst harmonic sum to be 

flM = ^r(0Ai(s) = - l)C(i! - 2j - 2 )k 2 „.(s - (2i + r) - 2). 

The Mellin transform of the second sum can be found in a completely analogous way: 
we have 

/2 (s) = (2"-^^ - l)C(s - 2j>2m(s - (2j + r)). 

Altogether, this yields the Mellin transform 


f*(s)=2f*is + 2)-f;is) 

s — 1. 




1)C(5 - 2j)K2m{s - {2j + r)). 


As in the proof of Lemma [3^ the growth conditions necessary for application of the 
converse mapping theorem il Theorem 4] hold. 

In order to analyze the poles of f*{s), we need to distinguish three cases, as 
C(s — 2j) has a simple pole at s = 2j + 1 and K 2 m(s — {2j + r)) has a simple pole 
at s = 2j + r + 1. The poles of r(s/2) at even s < 0 are canceled by the zeros of 
C(s —2j), unless s = j = 0. In that case, the pole is canceled by the factor (2^“^-^ — 1). 

First, let r = j = 0. Then, f*{s) has a simple pole at s = 1, because one of the 
poles of C(s) or K 2 m(s) cancels against the zero of (s — 1). Here, the residue of f*{s) 
is given by which translates to a contribution of y/rmr/A. This proves (18). 

Second, for r = 0 and j > 0, the function f*{s) has a pole of order 2 at s = 2j + 1. 
By expanding all the occurring functions, we hnd the Laurent expansion 


r{s) 



1 

(.-(2j+l))2 

1 

(s-(2j+l))2 


+ 


+ 


5V’(i+|)+27+log2+T 


s-(2j+l) 

ip(i+i)+27+31og2-2+An 


s-(2j+l) 


^ + 0 ( 1 ) 

+ 0 ( 1 ) 


for m G N, 
for m ^ N, 
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where 'ip{s) is the digamma function (cf. [5l 5.2.2|, see [5l §5.4(ii)| for special values). 
As the pole of order 2 contributes the factor logm, and the pole of order 1 


gives (19) is proved. 


Finally, consider r > 0. In this case we have two separate single poles at s = 2j' +1 
and s = 2j + r + 1. Computing the residues gives the growth contribution 


jT {j + \) + (j + ^) r (j + ) ( 2" 




which proves (20) 


□ 


Fortunately, when explicitly computing the expansion, all the logarithmic terms 
cancel out and we obtain the same behavior for admissible paths of even and odd 
length. The following theorem summarizes our hndings. 

Theorem 3 (Asymptotic analysis of admissible random walks on Z). The probability 
that a random walk on TL is admissible has the asymptotic expansion 


_ 1 4 88 

^ “ n 3n2 ^ 45n3 SlSn^ ' 675^5 31185^6 


976 3488 276928 ^ / 1 

\n' 


The expected height of admissible random walks on Z is given by 

, , r- 3V^ 539^^ 50713^^ 1 

4 16^^ 5760^^ 4838407^ 


7/ ’ 


where \/2 'k‘^ jA ^ 1.96870, and the variance of Hn can be expressed as 


(23) vg„, 28C(3)-.^ ,^^ 224C(3)-9.^ 1792C(3) - 67.= 


48 


2880n 

107520C(3) - 41897r3 

120960n2 


+ 0 K , 




where (28(C(3) —7r^)/8 ~ 0.33141. Generally, the r-th moment is asymptotically given 
by 


( 24 ) 


A . ^ + 1 

IT\ 


(^)(2^+i - \)2-'^l\{r + 
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Moreover, if rj = h/y/n satisfies 6/\/log n < rj < ^/\ogn/2, we have the local limit 
theorem 


(25) 


Sh) 


^ ft) ^ ^ 2^ 

Qn A/n 


AV2 


'TTn 


^((2A; + 1)V - l)exp 


{2k + 1) V 


k>0 


(26) 


^ P I 2/^2 


Remark. Again, the two expressions for the limiting density y are equivalent, as can 
be seen by an application of the Poisson sum formula. 

Proof. Analogous to Theorem The asymptotic expansions were again computed 
with the help of SageMath [T5], and a corresponding worksheet (containing these 
computations as well as some numerical comparisons) can be found at http: // 
arxiv.org/src/1503.08790/anc/randoni-walk_ZZ. ipynb. □ 


Remark. As every simple symmetric random walk of length n on Z occurs with 
probability 2“*^, we know that the number of admissible random walks on Z is 2"’g„. 
Thus, an asymptotic expansion for the number of admissible random walks follows 
directly from (21) upon multiplication by 2". This is sequence A167510 in HU. 



Figure 3. Numerical approximation of Bn/2"^ 
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Furthermore, in the introduction we illustrated that admissible random walks are 
strongly related to bidirectional ballot sequences. Since every bidirectional ballot 
sequence of length n + 2 corresponds to an admissible random walk of length n on 
Z (i.e., Bn = 2^~‘^qn-2), we are able to prove Zhao’s conjecture that was mentioned 
in the introduction. 


Corollary 4.2 (Bidirectional ballot walks). The number of bidirectional ballot walks 
Bn of length n can be expressed asymptotically as 


+ 


10 


+ 


764 


+ 


4952 


63n4 4725n5 31185n6 


on 

+ 01 ^ 


n‘ 


( 27 ) 73„-2”(i- + ^ + 

In Figure]^ we compare the exact values of Bn/2"^ with the values obtained from 
the asymptotic expansion in (27). 
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